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A Hamiltonian eikonal model for the multiple production in high energy hadron 
hadron collisions is presented and worked out in the aim of providing a simple 
frame for many different observables of these processes. The Hamiltonian formu- 
lation ensures that it has the unitarity built in by construction and the eikonal 
approximation makes easier the discussion of the possible spacial inhomogcncitics 
of the hadrons. The experimental data which are examined are the inelastic cross 
section and the single and double inclusive cross sections. 



1 Description of the model 

A Hamiltonian model for the description of the multiple production process in 
hadron-hadron collision is presented with the main aim of bringing together 
different observables within a unique frame. A particular attention is given 
to those features of the inelastic processes that can give informations on the 
proton structure. 

The physical ingredients of this model are the following: 
A sharp distinction is made between the soft dynamics, which provides the 
binding of the partons in the hadrons and also the final hadronization of the 
shaken-off partons and the hard dynamics that causes the parton scattering. 
The hard collision gives a finite transverse momentum to the partons, which 
remains however small with respect to the typical longitudinal momentum. 
Hard rescattering is included, but not hard branching of the partons. 
The discrete quantum numbers spin and colour are not taken into account. 

The model describes the hadrons as sets of bound partons which, due to the 
interaction, may become finally unbound and show up as jets; the hadroniza- 
tion process is not described. 

The only detailed kinematics is the transverse one, the whole treatment 
lies within the frame of the eikonal formalism. [1,2, 3] The sharp distinction 
between backward and forward degrees of freedom allows the following for- 
mulation: There are operators for the bound backward and forward partons 
a,b,af and operators for the unbound partons q,,c/. Both are local in b, 
the transverse impact parameter of the parton, and they have the standard 
commutation relations, where every backward operator commutes with every 
forward operator and every a commutes with every c; so we can write a free 
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Hamiltonian: 

H ° = E u I d 2 b[at(b)a v (b) + 4(b)c(b)] (1) 

v=f,b J 

The interaction that we want to describe is the hard collision of two bound 
partons that give rise to two unbound partons is such a way however that they 
keep their property of being either backward or forward. Thus the interaction 
Hamiltonian is written as: 

Hi = \Jd 2 bh b (b)h f {h) h v (b) =4(bK(b) + <4(b)c„(b) ■ ( 2 ) 

With this choice it results [H ,H/] = , but the theory is not trivial even 
though it has been much simplified, the S-matrix can be written in the form 
S = exp[— iHt] where r is an interaction time. 

An alternative, and more realistic, form gives a finite size to the hard 
interaction and a discretization of the transverse plane. The size A is related to 
the cut-off in the transverse momentum that must be put in order to be allowed 
to perform perturbative calculations, so the natural choice is A w l/p\; this 
choice lead also to the interpretation of r s=s l/p± w V~A. The commutation 
relations become [A v j, J = Sij6 UiV , and so on. In this way we get: 

n = J2ml, j A v , j + cl j c VJ ] (3) 

v,j 

Hi - (g/VA) ]T H bj ■ H fd H vJ = C^A^ + A\.C V ^ . (4) 
j 

The coupling constant g is dimensionless, it is related to the previous coupling 
constant by g = A/\/A. Correspondingly the S-matrix is 

S-Jp, , Sj=exp[-i(g/VZ)TH btj -H ftj ], (5) 

3 

with the previous interpretation of r, we get the simpler form: 

Sj = ex.p[-igH bJ ■ H ftj ] . (6) 

In order to apply this model one must choose a definite initial state; it will 
be factorized in the same way as the S-matrix: as far as its structure in a 
site j is concerned there are no strong indications. A possible choice, related 
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to some theoretical ideas about the nonpcrturbative partonic structure of the 
hadrons,[4,5] is a local coherent state, so we write 

\I >= II l J >i > l J >i= exp[-(|F 6 | 2 +|J>| 2 )/2] eMF b Al+F f A\] \ > j (7) 

3 

It has to be noted that the weight F of the coherent state may vary from site 
to site. For simplicity the index j will be not written out, whenever possible. 

It is easer to express \I >j and especially to perform the subsequent cal- 
culations in the basis generated by the auxiliary operators P and Q. 

P=(C + A)/V2 Q = (C - A)/V2 . (8) 

In terms of them on gets 

W = ^ ] ^[Py.j-PvJ + Qv,jQv,j] > H v ,j = -PyjPv.j ~ Qv.jQv.j (9) 

v,j 

2 Some results of the model 

2.1 Inelastic cross section 

In the discrete formulation for the states and for the S-matrix the inelastic 
cross section is now calculated. In the basis generated by the operators P and 
Q the operator Sj is diagonal, so it is easy to calculate the matrix element 
Sj = j < I\Sj\I >j, it has the expression: 

k\ ■••k4 

(10) 

The indices k\,k2,ks, refer to the quanta created byP^, Q\,pJ, Q| respec- 
tively. The normalizing factor is TV = cxp[— (\F b \ 2 + |F/| 2 )/2]. Using the 
representation: 

exp[— ig(ki— k2)(k^— k^)] = (2n)~ 1 J dudv exp[iuv+iau(ki— k2)+i(3v (ks— fct)] , 

(11) 

with a(i = g, the multiple sum in the expression of Sj can be transformed into 
an integral. In the final result we use the definitions T v = \F V \ 2 and we have: 

Sj — (27r) _1 / dudv exp[iuv] exp[— T(,(l — cos cm) — T/(l — cos{3v)] . (12) 
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When the distribution functions F v do not vary strongly from site to site 
one can devise a continuum limit. We consider a relation F w /\/A; therefore, 
if / is not singular, in the expression of Sj the terms \F\ 2 become small, the 
exponential in the integral representation can be expanded and integrated term 
by term with the result: 

Sjnl- \F b \ 2 \F f \\l - cosg) + (l/2)|F b | 2 |f>| 2 (|^| 2 + \F f \ 2 )(l cosg) 2 + . . . 

(13) 

With the normalization we are using the inelastic cross section at fixed hadronic 
impact parameter B is given by 

(7(B) = 2 < J|(l - 3tS)\I >-\< - S)\I >f . (14) 

The product of the matrix elements is expressed as the exponential of the sum 
of the logarithms and the sum A J2j i s finally converted into the integration 
J d 2 b with the final result: 

<7(B) = l-exp[-* J d 2 b\f b (b)\ 2 \f f (b-B)\ 2 +a 2 ---] (15) 

the parameter u — 2A(1 — cosg) has the role of elementary partonic cross 
section. 

The form of the inelastic cross section is quite usual, it contains however 
the explicit indication of the possible corrections due to rescattering, they will 
be discussed below, where a nonuniform model of the hadron will be explored. 
The cross section arises from the integration over the impact parameter, the 
result depends very much on the properties of the distribution functions /, 
which appear through their squared absolute values t v (b) = \f v (b)\ 2 , giving 
the transverse density of bound partons. In a first discussion the distribution 
is taken to be completely uniform in b by setting: 

t b (b)=p b #(R-\b\) , t f (b) = Pf 0(R-\b\); (16) 

elementary geometrical considerations give \B\ = 2R cos 7/2 with < 7 < n. 
The exponent in the integrand is given by the partial superposition of the 
two disks. Since the superposition area is W — R 2 {"f — sin 7) = nR 2 ^ The 
expression for the inelastic cross section is: 

PIT 

a in — 2-kR 2 I d7sin7 1 — exp[— v£] , (17) 

Jo 

with v = apbPfTrR 2 . 
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It is possible to give a simple analytical form for a in in the two limiting 
situations of very small or very large v. In the first case one gets 



= ttR 2 ■ v . 



(18) 



In the second case we can start from the expression 



2ttR 2 [2 - D{v)} 



(19) 



where the real function D(v) defined by eq(17) is monotonically decreasing, 
for large v it results D(v) w 2(&v 2 /n 2 )^ 1 ^ ■ T(2/3) so that the geometrical 
limit of black disks AitR 2 is approached. 

2.2 Pair and double-pair production 

The production of a pair can seen either as production of a backward parton or 
of a forward parton, the Hamiltonian being fully symmetric, so we can choose, 
arbitrarily to look at the forward particles; successively we shall investigate 
how much the rescattering processes may destroy the sharp correlation between 
backward and forward scattered partons. We start from the computation of 
the inclusive production from a single site. Straightforward, although lengthy 
calculations, easier in the basis generated by the operators P and Q, yield the 
result 



< Xj >= j < I\Sjc}C f Sj\I >j= (T f /2) [1 - exp[-T 6 (l - cos(2 5 )]] . (20) 



We can now go to the continuum limit, always under the hypothesis of 
smooth distributions t v (h); we find in eq(20) the function cos(2g) instead of 
cos g, so we must define a quantity n— (A/2)[(l — cos(2g)], which is related to 
a in this way: k = a ■ [1 + cos(g)]/2. The two constants coincides evidently for 
small g, when both are: a = k = g 2 A, but the rescattering corrections, higher 
powers in g 2 , are different. 

In the smooth continuous limit we expand the exponential of eq. (4.3) and 
we get the usual expression: 




(21) 



and the inclusive one-particle cross section is 




(22) 
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Analogously the two-particle inclusive cross section is found to be 



D 2 = k 2 J d 2 Bd 2 bd 2 b't f {b)t b {b - B)tf(b')t b (b' - B); . (23) 

A ratio of the quantities that have been now calculated which is of phe- 
nomenological interest is a e g = [Di] 2 /D 2 ,[6,7] In the limit of rigid disk, using 
the geometrical considerations one gets 

°*-T=ms?) (24) 

The previously calculated expression of a in is really the hard part of the 
total inelastic cross section, where as "hard" part we intend the contribution 
of all the events where at least one hard scattering happens. If we believe 
that, in going on with the total energy these events become more and more 
important we would like to have this term not too small with respect to the 
experimental <7j„, which in turn appears to be sizably larger, of about a factor 
2, with respect to a e g, so in this model we would like to approach, even though 
not reach, the black- hadron, limit which produces o in — AirR 2 . 



2.3 Multiplicity distribution and forward-backward correlations 



The distribution of the multiplicities of the produced pairs is calculated by 
defining the projection operator over the number of free partons. Since we 
have always the sharp distinction between backward and forward particles we 
can choose one of the two particle to define the produced pair: to be definite 
we take the forward particle as signal of the pair production. For a fixed site 
j the number projector may be expressed as: 

V n = -.: [C^C] n e' c ' ic : (25) 
n! 

The colon indicates the normal ordering of the C-operators. From this defi- 
nition one can perform calculations whose qualitative results can be stated in 
this form: when the total production is not very copious the result may be 
written as the sum of two Poissonian distributions, one reflects the incoming 
coherent state the other the rescattering effect , if however the production in 
a single site is high so that the rescatterig is very important the resulting ex- 
pression may be put into the form of a Poissonian distribution times another 
factor, but this further factor is not a small correction, it changes in essential 
way the shape of the distribution. 
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The same qualitative result is obtained by studying the forward-backward 
correlations; we can calculate both the variance £/ =< X 2 > — < Xf > 2 with 

< Xj >=< I\S^ C)-C f C]C f S\I > and the covariance 

S /j6 =< V > - < X f >< X b > with: < V > = < I\S^C)C f ClC b S\I > 
The actual calculations show, as expected, that the correlation 



for small values of T v , and so for small production, goes to 1, but when T v 
becomes very large it goes to zero. 

2.4 Non uniform hadrons 

We wish now to explore the possibility that the hadron, and its projection 
over the transverse plane, shows strong inhomogcncitics in the matter density. 
This is represented by assuming that there are black spots, that cover a limited 
amount of the transverse area, while a much fainter "gray" background fills 
uniformly the rest of the hadron. 

\h >j= [zcxp[-|F 6 | 2 /2] exp[F b 4] + y exp[-|G 6 | 2 /2] cxp[G 6 4]] | > (27) 

Since the two coherent states are not orthogonal the normalization condition 
is complicated, however if the two thickness are very different the cross term 
in normalization condition is very small and we are left with \x\ 2 + \y\ 2 ~ 1 . 
The expression for the S-matrix may be given in the form 

S = S S C , S = S W , w = W/A (28) 

the factor S W gives the contribution of the background scattering. 

When one uses the assumption that y « x i.e. that the spots cover a 
small part of the transverse area it is possible to give a simple expression to 
S c 

S c = 1 - wy 2 A[ Pb (l - a Pf ) + p f (l - a Pb )] (29) 

It may be noticed that the factor, wy 2 A, represents the part covered by black 
spots within the interacting area of the hadrons at given B. 

The basic elements for calculating the pair and double pair production have 
already been given. The local production amplitude is the sum of 4 terms, it 
will be indicated as: < Xj >= [x i X OM + x 2 y 2 (X s , + X a , s ) +y A X s ^ s ] . The first 
term represents the pure background interaction, the second and third terms 
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represent the spot background interaction, the fourth term gives the spot-spot 
interaction. When one consider explicitly the superposition of the two disks it 
appears that the four terms have the same geometrical properties and we get 
for er c ff the same expression as in the case of a uniform hadron. So, at first 
sight it seems that nothing is gained by introducing an inhomogeneity into the 
hadron, but in fact some new features are present. The expression of a c g is 
purely geometrical, it does not contain the parameters x,y. The expression 
of a in can be obtained from the S-matrix (see eq.28, 29), so it depends on y, 
more in general both on x and on y. 

In order to give a clear, although unrealistic evidence of this fact we can 
consider the limit in which the background is so thin that it contribute neg- 
ligibly in the inelastic cross section, then the S-matrix element, depends only 
on the spot-spot interaction In this situation the inelastic cross section goes 
to zero with y 4 , the corresponding picture would be: few, wholly black spots 
distributed in a wide and very thin background. The conclusion of this analy- 
sis is therefore that in order to have a c s < Oi n the hadron should be compact, 
i.e. without holes or transparent regions and also with quite sharp edges. 

3 Further considerations and conclusions 

The model allows further derivations that are here only mentioned: it is pos- 
sible to introduce correlations in b in the case of nonuniform hadron; it is 
also possible to treat to some extent the longitudinal degrees of freedom pro- 
vided the distinction between forward and backward particles is kept valid. 
Moreover, in some simple cases analytically and in more general situation in 
numerical way, it is possible to go beyond the two extreme cases that have been 
illustated here i.e. the weakly interacting or the totally absorbing situation. 

The model presented allows a systematization of different aspects of the 
hard processes in multiparticle production with a particular attention to the 
unitarity corrections and suggests also some interpretations in terms of hadron 
structure. The connection with QCD is not direct as it appears from the fact 
that the interaction term is quartic while the fundamental QCD interaction 
term is cubic, in other words the fundamental input is the parton hard scat- 
tering, not the branching process. 
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